arXiv: 1503.05983vl [cond-mat.mtrl-sci] 20 Mar 2015 


GEOPHYSICAL RESEARCH LETTERS, VOL. ???, XXXX, DOI:10.1029/, 


Thermoelasticity of Fe^+-bearing bridgmanite 

Gaurav Shukla/ Zhongqing Wu,^ Han Hsu,^ Andrea Floris,"^ Matteo 
Cococcioni,^ Renata M. Wentzcovitch/’^ 


Corresponding anthor: Ganrav Shnkla, School of Physics and Astronomy, University of Min¬ 
nesota, Minneapolis, Minnesota, USA. (shukla@physics.umn.edn) 

^School of Physics and Astronomy, 

University of Minnesota, Minneapolis, 

Minnesota, USA 

^School of Earth and Space Sciences, 

University of Science and Technology of 
China, Hefei, Anhui, China 

^Department of Physics, National Central 
University, Jhongli District, Taoyuan 32001, 

Taiwan 

'^Department of Physics, King’s College 
London, England, United Kingdom 

^Theory and Simulation of Materials 
(THEOS), Ecole polytechnique federale de 


DRAFT 


March 23, 2015, 12:35am 


DRAFT 



X - 2 


SHUKLA ET AL.: THERMOELASTICITY OF IRON-BEARING BRIDGMANITE 


We present LDA+U calculations of high temperature elastic properties of 
bridgmanite with composition (Mg(i_ 3 ;)Fe^+)Si 03 for 0 < x < 0.125. Re¬ 
sults of elastic moduli and acoustic velocities for the Mg-end member (x=0) 
agree very well with the latest high pressure and high temperature exper¬ 
imental measurements. In the iron-bearing system, we focus particularly on 
the change in thermoelastic parameters across the state change that occurs 
in ferrous iron above ~30 GPa, often attributed to a high-spin (HS) to in¬ 
termediate spin (IS) crossover but explained by hrst principles calculations 
as a lateral displacement of substitutional iron in the perovskite cage. We 
show that the measured effect of this change on the equation of state of this 
system can be explained by the lateral displacement of substitutional iron, 
not by the HS to IS crossover. The calculated elastic properties of (Mgo.875Feo)i25)Si03 
along an adiabatic mantle geotherm, somewhat overestimate longitudinal ve¬ 
locities but produce densities and shear velocities quite consistent with Pre¬ 
liminary Reference Earth Model data throughout most of the lower mantle. 
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1. Introduction 

Thermoelastic properties of lower mantle minerals provide a direct link to seismolog- 
ical observations. In order to constrain the composition and thermal structure of the 
Earth’s lower mantle, a detailed investigation of elastic properties of the constituent 
minerals is needed. Bridgmanite, (Mg,Fe,Al)(Si,Fe,Al )03 perovskite (Pv), is the main 
constituent of the lower mantle along with ferropericlase, (Mg,Fe)0, CaSiOs perovskite, 
and (Mg,Fe,Al)(Si,Fe,Al )03 post-perovskite (PPv). Although there has been consider¬ 
able progress in measurements of elastic properties of these minerals at high pressures 
and temperatures [Yeganeh-Haeri, 1994; Sinnelnikov et al., 1998; Piquet et al., 2000; An- 
drault et al, 2001; Aizawa et al, 2004; Jackson et al, 2004, 2005b; Sinogeikin et al, 2004; 
Li and Zhang, 2005; Murakami et al., 2007; Lundin et al., 2008; Boffa Ballaran et al., 
2012; Chantel et al., 2012; Murakami et al., 2012; Dorfman et al., 2013], owing to ex¬ 
treme pressure and temperature conditions in the lower mantle the availability of data 
on Fe-bearing bridgmanite is quite limited to well constrain lower mantle composition 
and temperature. Several experiments and hrst-principles calculations have shown that 
pressure induced iron spin crossover [Badro et al, 2003, 2004; Tsuchiya et al, 2006; Lin 
et al, 2012] affects elastic properties of (MgFe)O ferropericlase (fp) [Goncharov et al, 
2006; Crowhurst et al., 2008; Marquardt et al, 2009; Wentzcovitch et al, 2009; Wu et al, 
2009; Antonangeli et al., 2011; Wu et al., 2013; Wu and Wentzcovitch, 2014]. In the case 
of iron-bearing bridgmanite, effects of spin crossover on elastic properties have been quite 
uncertain due to possible coexistence of ferrous (Fe^’*') and ferric (Fe^’*') iron along with 
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the more complex perovskite crystal structure. Experimentally, it has been difficult to 
isolate the effects of ferrous and ferric iron unambiguously. 

Acoustic velocity measurements using an ultrasonic technique [Chantel et ai, 2012] and 
compression curve measurements [Boffa Ballaran et ai, 2012] in iron-bearing bridgman- 
ite, demonstrated a signihcant change in bulk modulus across the state change in iron. 
By htting experimental data with S^Aorder Birch-Murnaghan equation of state to small 
(0-40 GPa) and large (0-75 GPa) pressure ranges, Boffa Ballaran et al. [2012] found ambi¬ 
ent condition bulk modulus smaller (245 GPa) and larger (253 GPa), respectively. Since 
these htting pressure ranges had no effect on the iron-free phase, Boffa Ballaran et al. 
[2012] attributed this behavior to change in compression mechanism caused by the high- 
spin (HS) to intermediate-spin (IS) crossover [McCammon et ai, 2008]. However the HS 
to IS crossover has not been found in hrst-principle studies [Bengtson et ai, 2009; Hsu 
et ai, 2010, 2011; Hsu and Wentzcovitch, 2014]. Using the local density approximation 
augmented by the Hubbard type correction (LDA-I-U) method, Hsu et al. [2010] showed 
that iron in Fe^’''-bearing bridgmanite should remain in the HS state throughout the lower 
mantle and should, instead, undergo a pressure induced displacement producing a state 
with increased iron Mdssbauer quadrupole splitting (QS). The calculated low (1.9-2.4 
mm/s) and high (3.3-3.9 mm/s) QS states were consistent with experimental measure¬ 
ments [Lin et ai, 2012; McCammon et al., 2008, 2013; Potapkin et ai, 2013; Kupenko 
et ai, 2014]. Owing to lack of sufficient experimental measurements on elasticity of iron¬ 
bearing bridgmanite, effects of the proposed HS to IS crossover [McCammon et ai, 2008], 
or iron displacement [Bengtson et al, 2009; Hsu et ai, 2010, 2011; Hsu and Wentzcovitch, 
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2014], have not been properly understood. To clarify this issue we have calculated the high 
temperature and high-pressure elastic moduli and acoustic velocities of (Mgi_ 3 -Fe 3 ;)Si 03 
for x=0 and x=0.125, with iron in low- and high-QS states. We also compare the effects 
of the HS to IS transition [Hsu and Wentzcovitch, 2014] in the static equation of state of 
Fe^’''-bearing bridgmanite. 

Unlike in (MgFe^+)(SiFe ^+)03 [Hsu et ai, 2011], the high-spin (HS) to low-spin (LS) 
transition in (MgFe^^)Si 03 has not been found in the lower-mantle pressure range using 
the LDA-I-U method [Hsu et ai, 2010; Hsu and Wentzcovitch, 2014]. However, there 
have been several studies suggesting that HS to LS transition might occur at around 
~ 90 GPa [Badro et al, 2004; Li et al, 2004, 2006; Jackson et al, 2005a; McCammon 
et al, 2008; Umemoto et al., 2008; Lin et al., 2012] and might have consequences for the 
vibrational spectrum [Caracas et al, 2014]. Recently Zhang et al. [2014] found that iron¬ 
bearing bridgmanite dissociates into nearly iron-free phase and into an iron-rich hexagonal 
silicate phase in the deep lower mantle. As the dissociation phase boundary is still not 
well known, we present results for elasticity of Fe^^-bearing bridgmanite with constant 
iron concentration in the entire pressure range of the lower mantle. 

2. Method and Calculation Details 

We have used density functional theory (DFT) within the local density approximation 
(LDA) [Ceperley and Alder, 1980] augmented by the Hubbard U (LDA-I-U) calculated 
self-consistently [Cococcioni and de Gironcoli, 2005; Kulik et al, 2006; Hsu et al, 2009; 
Leiria Campo Jr. and Cococcioni, 2010; Himmetoglu et al, 2014]. Ultrasoft pseudo¬ 
potentials [Vanderbilt, 1990] have been used for Fe, Si, and O. For Mg, a norm-conserving 
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pseudo-potential generated by von Barth-Car’s method has been used. A detailed descrip¬ 
tion about pseudo-potentials is available in Umemoto et al. [2008]. The plane-wave kinetic 
energy and charge density cut-off are 40 Ry and 160 Ry, respectively. Structural optimiza¬ 
tion at arbitrary volumes has been performed using variable cell-shape damped molecular 
dynamics [Wentzcovitch, 1991; Wentzcovitch et al, 1993]. Self- and structurally-consistent 
Usc=3.1 eV, reported by Hsu et al. [2010] using the linear response approach [Cococcioni 
and de Gironcoli, 2005; Kulik et al., 2006], has been used for high-spin low-QS and high- 
QS iron. Vibrational density of states (VDoS), needed to calculate free energy within 
quasi-harmonic approximation (QHA) [Carrier et al., 2007], have been calculated using 
density functional perturbation theory (DFPT) [Baroni et al., 2001] using the LDA-I-U 
functional [Floris et al., 2011]. (Mgi_ 2 ,Fea;)Si 03 for x=0 and 0.125 have been investigated 
in a 40-atom super-cell. Elastic properties for 0<x<0.125 are linearly interpolated using 
x=0 and x=0.125 results. Electronic states have been sampled on 2 x 2 x 2 Monkhorst- 
Pack k-point grid shifted by (1/2,1/2,1/2) from the Brillouin Zone origin. VDoS have 
been obtained by calculating dynamical matrices on a 2 x 2 x 2 q-point grid and inter¬ 
polating the obtained force constants on a 8 x 8 x 8 q-point grid. The structure for x=0, 
and x=0.125 with high-QS state have been optimized at twelve pressures between -10 
and 150 GPa. The low-QS structure has been optimized up to 60 GPa only. Beyond this 
pressure unstable phonons appear. To obtain static elastic coefficients (zero temperature), 
Cij, at each pressure, positive and negative strains of 1% magnitude have been applied 
and after relaxing the internal degrees of freedom, associated stresses have been calcu¬ 
lated. Thermoelastic moduli have been calculated using a semi-analytical approach [Wu 
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and Wentzcovitch, 2011]. This method uses an analytical expression for strain Griineisen 
parameters to calculate the thermal contribution to elastic coefficients using the quasi¬ 
harmonic approximation. This method allows calculations of thermal elastic coefficients 
using static values and vibrational density of states for unstrained conhgurations only. It 
is almost two orders of magnitude more efficient than the fully numerical method [Karki 
et ai, 1999; Wentzcovitch et al., 2004], which required vibrational density of states also 
for strained conhgurations. It also appears to be more accurate since it avoids numerical 
differentiation on a reduced number of grid points as well as calculations of VDoS and free 
energies for strained conhgurations. The method has been applied successfully to several 
minerals already [Wu and Wentzcovitch, 2011; Nunez-Valdez et al, 2012a, b, 2013; Yang 
and Wu, 2014]. Calculations have been performed using VLab cyberinfrastructure at the 
Minnesota Supercomputing Institute [da Silveira et al, 2008]. 

3. Results and Discussion 

3.1. Elasticity of Iron-free and Fe^+-bearing Bridgmanite 

The calculated aggregate elastic moduli, K and G, and acoustic velocities, Vp and Vs, of 
iron-free bridgmanite compare well with Brillouin scattering measurements of Murakami 
et al. [2007, 2012] as shown in Fig. 1(a) and 1(b). Within experimental uncertainties, 
ultrasonic measurements of Li and Zhang [2005] and Chantel et al. [2012] also compare 
well with our results. Table 1 compares calculated elastic properties at zero pressure 
and 300 K with measurements and simulations at ambient conditions. As compared to 
experimental measurements, our calculated volume is overestimated by ~ 0.5% while 
elastic moduli and acoustic velocities are underestimated by ~ 1-3% , which is typical for 
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LDA calculations after including zero-point energy [Karki et al, 1999, 2001; Wentzcovitch 
et al, 2004, 2010a, b]. Calculated pressure derivatives of elastic moduli, K' and G', 
are 3.96 and 1.79, respectively. Our K' agrees reasonably well with previous studies 
while G' compares well with those reported by Sinnelnikov et al. [1998] and Zhang et al. 
[2013]. Comparing acoustic velocities and elastic moduli at higher temperatures, especially 
shear velocities and shear moduli at 2,700 K by Murakami et al. [2012], we find that the 
temperature dependence of elastic properties has been captured well. As shown in Fig. 
2, our calculated and ^ are smaller in magnitude compared to values obtained by 
ultrasonic [Sinnelnikov et al, 1998; Aizawa et al, 2004; Li and Zhang, 2005] at 0 GPa and 
Brillouin scattering measurements [Chantel et al., 2012] at 22 GPa. Calculated values of 
compare well with those of previous calculations [Oganov et al, 2001b; Marion and 
Cohen, 2002; Wentzcovitch et al, 2004; Zhang et al, 2013]. Calculated ^ agrees well 
with Brillouin scattering measurements by Murakami et al. [2012] and non-self consistent 
calculation by Marion and Cohen [2002] but differs substantially from that of Oganov 
et al. [2001b], Wentzcovitch et al. [2004], and Zhang et al [2013]. 

The semi-analytical method used here had produced well the temperature dependence 
of elasticity for all applied materials [Wu and Wentzcovitch, 2011; Nunez-Valdez et al, 
2012a, b, 2013; Yang and Wu, 2014]. This method has also produced slightly better 
results for MgO [Wu and Wentzcovitch, 2011] than the fully numerical method [Karki 
et al, 1999; Wentzcovitch et al, 2004]. This is because this method avoids calculations 
of vibrational frequencies in strained configurations and numerical errors associated with 
hnite difference calculations of elastic coefficients. 
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The elastic moduli and acoustic velocities of iron-bearing bridgmanite at lower mantle 
pressures are also shown in Fig. 1. To compare with measurements of Boffa Ballaran 
et al. [2012] and Chantel et al. [2012] in the iron-bearing phase, we estimated elastic 
properties for 0<x<0.125 by linearly interpolating our calculated results for x=0 and 
x=0.125. The calculated volume at ambient conditions is in agreement with measurements 
by Boffa Ballaran et al. [2012] for x=0.04 and is smaller by ~2% when compared with 
measurements by Chantel et al. [2012] for x=0.05 (Table 1). As shown in Fig. 1(b), 
comparison of calculated acoustic velocities, Vp and Vs, for x=0.05 with the only available 
experimental measurements for (Mgo, 95 Feoj) 4 Feoj)i)Si 03 [Chantel et al., 2012] is also good 
except for the shear velocity data in the lower pressure range. As seen in Fig. 1 (c) and 
(d) and in Table 1, the presence of iron has very little to no effect on bulk modulus, K5, 
but shear modulus G softens by ~2 to 6%, consistent with previous static hrst-principles 
studies [Kiefer et al, 2002; Lundin et al, 2008]. The effect of iron on the temperature 
dependence of acoustic velocities and elastic moduli is not very signihcant (see Fig. 1). 

3.2. Effect of Pressure Induced Iron Displacement on Elasticity 

We have investigated the elastic consequences of the pressure induced iron displacement 
[Hsu et al., 2010; Hsu and Wentzcovitch, 2014] by calculating thermal elastic moduli 
and acoustic velocities for x=0.125 with iron in low-QS and high-QS sites. The low-QS 
structure is stable up to 60 GPa in LDA-I-U calculations. Above this pressure unstable 
phonon modes start appearing. As shown in Fig. 3, vibrational density of states (VDoS) 
for x=0 and for x=0.125 with iron in the high-QS and low-QS sites are not very different 
except for a small shift towards lower frequencies owing to the presence of iron. However, 


DRAFT 


March 23, 2015, 12:35am 


DRAFT 


X - 10 


SHUKLA ET AL.: THERMOELASTICITY OF IRON-BEARING BRIDGMANITE 


there is an extra peak at ~ 150 cm“^ in the VDoS if iron is in the low-QS site. It is worth 
noting that the vibrational mode associated with this extra peak in the low-QS structure 
is infra-red active while the same mode is infra-red inactive in the high-QS structure. 
This information may provide a crucial test to investigate the state of ferrous iron in 
bridgmanite. 

The effect of iron displacement (low-QS to high-QS transition) on unit-cell volume, 
acoustic velocities, and elastic moduli of iron-bearing bridgmanite with x=0.125 is shown 
in Fig. 4. To quantify the magnitude of the property change, we dehne the contrast A in 
a particular property M as AM = ^^igh-Qs-Mio^.Qs ^ shown in Fig. 4, the volume 

^^^low—QS 

contrast is very small (~0.1%), consistent with Lundin et al. [2008] and Boffa Ballaran 
et al. [2012], and is unlikely to be detected in experiments. The contrast in bulk modulus. 
Kg, and acoustic velocities, Vp and Vs, are also small. However, the shear modulus 
contrast is about ~0.5-1.5%. This change in G, was also suggested by Bojfa Ballaran 
et al. [2012] using spontaneous strain analysis. To investigate the htting pressure-range 
behavior observed by Boffa Ballaran et al. [2012] and Chantel et al. [2012] the calculated 
results for low-QS and high-QS have been combined at 30 GPa, the pressure above which 
the HS to IS crossover was proposed by McCammon et al. [2008] and Potapkin et al. 
[2013]. The 3'’'^-order Birch-Murnaghan equation of state was used to ht these combined 
results. As shown in Table 2, depending on the pressure-range of the combined results, 
an isothermal bulk modulus softening of ~3 to 7 GPa has been found. This softening, 
however, reduces with increasing temperature. 
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Intermediate spin (IS) states of the iron in Fe^’''-bearing bridgmanite are energetically 
unfavored in the lower mantle pressure-range [Hsu et ai, 2010; Hsu and Wentzcovitch, 
2014] and phonon calculations have been extremely difficult for these states. To under¬ 
stand the overall trend of htting parameters of combined low-QS and IS data, we have 
used static data for low-QS state and IS state with QS value 0.9-1.6 mm/s [Hsu and 
Wentzcovitch, 2014]. The static bulk modulus at zero pressure for low-QS and IS states 
are 263.2 GPa and 254.1 GPa, respectively (Figure S3). When static low-QS and IS re¬ 
sults have been combined at 30 GPa and htted within different pressure-ranges of 0-150 
GPa and 0-120 GPa, static bulk modulus of the combined ht are 305 GPa and 348.7 
GPa, respectively (Figure S4). The htting pressure-range trend of combined low-QS and 
high-QS is the same as that observed by Boffa Ballaran et al. [2012] and Chantel et al. 
[2012] while of combined low-QS and IS states the htting trend is opposite. These hndings 
substantiate our understanding that ferrous iron in perovskite remains in the HS state in 
the lower mantle pressure range, but instead, undergoes a horizontal displacement in the 
perovskite cage which results in an increase in Mdssbauer quadrupole splitting. 

4. Geophysical Significance 

As shown in Fig. 5, elastic moduli {Ks, G), acoustic velocities {Vp, Vs), and densi¬ 
ties (p) are calculated along the Boehler’s geotherm [Boehler, 2000] using new [Wu and 
Wentzcovitch, 2011] and previous [Karki et al, 1999; Wentzcovitch et al, 2004] methods 
and the results are compared with seismic values from the Preliminary Reference Earth 
Model (PREM) [Dziewonski and Anderson, 1981]. Owing to the smaller magnitude of ^ 
obtained here, present results for G and Vs diher more from previous results by Wentz- 
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covitch et al. [2004] in the deep lower mantle region, where temperatures are higher. This 
difference is sufficient to alter conclusions regarding mantle composition. Inclusion of 
12.5% Fe^'*' in MgSiOs perovskite produces good agreement of p, V 5 , and G with PREM 
values (Fig. 5c,d). Therefore, comparison of calculated Vs with PREM’s Vs may indeed 
suggest a more perovskitic lower manle [Murakami et al, 2012]. However, relatively large 
values of Ks and Vp are found, suggesting that the lower mantle may accommodate a 
reasonable amount of ferropericlase, (MgFe)O, and CaSiOs pereovskite as previously in¬ 
dicated [Wentzcovitch et al, 2004]. Resolving whether the lower mantle composition is 
pyrolytic or perovskitic still needs more detailed investigations of the effect of AI 2 O 3 and 
Fe 203 on the elasticity of bridgmanite along with the elasticity of CaSi 03 perovskite. Ul¬ 
timately, the analysis of one dimensional velocity prohles might be too limited to address 
this question more conclusively. Analysis of lateral heterogeneities, including effects of 
spin crossovers in (MgFe)O [Wu and Wentzcovitch, 2014] and in bridgmanite might prove 
more useful to address this question. 

5. Conclusions 

We have presented hrst-principles LDA-I-U calculations of aggregate elastic moduli and 
acoustic velocities for Fe^'''-bearing bridgmanite at lower mantle conditions. Using the new 
semi-analytical approach for high temperature elasticity [Wu and Wentzcovitch, 2011], 
we hud an unexpected difference of temperature gradients with respect to previous fully 
numerical calculations [Wentzcovitch et al., 2004] also based on the quasi-harmonic ap¬ 
proximation. New results for ^ and for iron-free bridgmanite are in good agreement 
with the latest high temperature high pressure measurements of these quantities by Mu- 
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rakami et al. [2012], which also differ from previous calculations and previous experiments. 
Overall, it appears that pressure and temperature gradients of bulk and shear moduli are 
better estimated now. Calculated acoustic velocities for iron-bearing bridgmanite also 
agree well with ultrasonic measurements of Chantel et al. [2012]. We do not find any 
significant effect of iron-site change (low-QS to high-QS) on sound velocities and elastic 
moduli, except for ~1% increase in shear modulus, consistent with findings by Boffa Bal- 
laran et al. [2012]. Effects of fitting combined data for low-QS and high-QS states on 
elastic properties are consistent with previously observed fitting pressure-range behavior 
observed by Boffa Ballaran et al. [2012] and Chantel et al. [2012]. Calculated shear mod¬ 
ulus and shear velocities of bridgmanite along a typical mantle geotherm are qualitatively 
more consistent with PREM values now compared with previous first principles results 
[Wentzcovitch et al, 2004]. Differences in bulk modulus and longitudinal velocities remain 
noticeable. 
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Figure 1. (Color online) Pressure and temperature dependence of elastic moduli and sound 
velocities for (Mgi_ 3 ;Fe^''')Si 03 with x=0 and 0.05 (a, b), and of dM/dx (M= K 5 , G, Vp, Vs) 
(c, d). Solid (dashed) lines represent first-principles results within (outside) the validity of quasi¬ 


harmonic approximation. BS: Brillouin scattering, US: Ultrasonic technique. 


DRAFT 


March 23, 2015, 12:35am 


DRAFT 
















SHUKLA ET AL.: THERMOELASTICITY OF IRON-BEARING BRIDGMANITE 


X- 25 



T3 

O 


0 


4000K 


(a) 


Simulations at T|yy|=2500K 
Oganov et al. [2001b] 
Marton and Cohen [2002] 
Wentzcovitch et al. [2004] 
Zhang et al. [2013] (LDA) 
Zhang et al. [2013] (GGA) 



(b) 


V 

o 


Experiments 
Sinnelnikov et al. [1998] at 300K 
Aizawa et al. [2004] at 300K 
Li and Zhang [2005] at 300K 
Chantel et al. [2012] at T|y|=500K 
Murakami et al. [2012] at 7^=1500K 


0 30 60 90 

Pressure (GPa) 


120 


150 


Figure 2. (Color online) Pressure and temperature dependence of (a) and (b) ^ for iron- 
free bridgmanite. Results from this study (lines) are compared with previous simulations (black 
symbols) and experimental data (red symbols). The values of and ^ obtained from previous 
simulations [Oganov et ai, 2001b; Marton and Cohen, 2002; Wentzcovitch et al, 2004; Zhang 
et ai, 2013] and experimental measurements [Chantel et ai, 2012; Murakami et ai, 2012] shown 
here are calculated at the midpoint between extreme temperatures [i.e., ^ 

where Tm = and M=(K 5 , G)]. Error bars on Oganov et al. [2001b] and Marton and Cohen 
[2002] results are reproduced from Zhang et al. [2013] while on Wentzcovitch et al. [2004] and 
Zhang et al. [2013] they are smaller than the symbols size. 
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Figure 3. (Color online) Vibrational density of states (VDoS) at P = 0 for iron-free bridg- 
nianite using LDA, and for low and high QS state of (Mgo.875Feo)i25)Si03 using LDA-I-U. 
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Figure 4. (Color online) Contrast (AM) in (a) unit-cell volume, (b) compressional and shear 
velocity (Vp, 14')) and (c) bulk and shear modulus {Ks, G) between low-QS and high-QS states 
in (Mgo.875Feo:i25)Si03. 
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x=0 x=0.125 



Figure 5. (Color online) Elastic moduli {Ks, G), sound velocities {Vp, 1 / 5 ), and density (p) for 
(Mgi_ 3 ;Fe^’'')Si 03 with x=0 (a, b) and x=0.125 (c, d) calculated along the Boehler’s geotherm 
[Boehler, 2000] using present [Wu and Wentzcovitch, 2011] and previous [Wentzcovitch et al, 
2004] methods. 
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Table 1. Volume, sound velocities, elastic moduli and their pressure derivatives at ambient 


pressure and temperature conditions. BS: Brillouin scattering, US: Ultrasonic, XRD: X-ray 


diffraction. 



X 

Vp{km/s) 

Vs{kmj a) 

As(GPa) 

G(GPa) 

K' 

G' 

References 

163.2 

0 

10.72 

6.42 

245.5 

168.2 

3.96 

1.79 

This study, single crystal/LDA 

163.3 

0.05 

10.63 

6.34 

245.8 

167.02 

4.03 

1.79 

This study, single crystal/LDA-|-U (Low- 
QS) 

162.3(0.6) 

0 

11.04 

6.57 

264(5) 

177(4) 

- 

- 

Yeganeh-Haeri [1994], sinlge-crystal/BS 

162.3 

0 

- 

6.51 

- 

176(5) 

- 

1.8(0.4) 

Sinnelnikov et al, [1998], poly-crystal, US 

162.4(0.5) 

0 

10.84(0.1) 

6.47(0.06) 

253(5) 

173(3) 

- 

- 

Sinogeikin et al. [2004], poly-crystal/BS 

162.4(0.5) 

0 

10.88(0.06) 

6.54(0.03) 

253(3) 

175(2) 

- 

- 

Sinogeikin et al. [2004], sinlge-crystal/BS 

162.3 

0 

10.86(0.05) 

6.49(0.04) 

253(2) 

173(1) 

4.4(0.1) 

2.0(0.!) 

Li and Zhang [2005], poly-crystal/US 

- 

0 

10.85(0.03) 

6.49(0.03) 

- 

172.9(1.5) 

- 

1.56(0.04) 

Murakami et al. [2007], poly-crystal/BS 

162.2(0.2) 

0 

10.82(0.06) 

6.54(0.04) 

247(4) 

176(2) 

4.5(0.2) 

1.6(0.1) 

Chantel et al. [2012], poly-crystal/US 

166.5(0.2) 

0.05 

10.60(0.06) 

6.46(0.04) 

236(2) 

174(1) 

4.7(0.1) 

1.56(0.5) 

Chantel et al. [2012], poly-crystal/US 

162.4 

0 

- 

- 

267 

180 

4.10 

- 

Oganov et al. [2001b], Simulations 

162.49 

0 

- 

- 

269(2) 

159(2) 

- 

- 

Marton and Cohen [2002], Simulations 

162.4 

0 

- 

- 

250.5 

172.9 

4.01 

1.74 

Zhang et al. [2013], Simulations 

Compression 

Studies 

v{A^) 




XT’ 


K' 



162.3 

0 



259.5 


3.69 


Piquet et al. [2000] 

162.4 

0 



248.8 


4 


Andrault et al. [2001] 

162.7 

0.05 



255.4 


4 


Andrault et al. [2001] 

164.1 

0 



247 


3.97 


Karki et al. [2001], Simulations 

162.30 

0 



255-261 


4 


Lundin et al. [2008] 

162.18 

0.09 



257-259 


4 


Lundin et al. [2008] 

163.30 

0.15 



257-259 


4 


Lundin et al. [2008] 

163.09(0.06) 

0.04 



253(2) 


3.99(0.07) 


Boffa Ballaran et al. [2012], (0-75 GPa) 

163.16(0.06) 

0.04 



245(4) 


4.4(0.03) 


Boffa Ballaran et al. [2012], (0-40 GPa) 

166.5(0.2) 

0.05 



246(2) 


4 


Chantel et al. [2012], 

163.0 

0.09 



251(13) 


4 


Dorfman et al. [2013] 


Table 2. Calculated equation of state parameters for Low-QS, High-QS, and combined results 
at 30 GPa for low-QS and high-QS fitted with different pressure range. 



300K 

700K 


v{A^) 

Kt 

K' 

U(i3) 

Kt 

K' 

Low-QS 

163.81 

247.16 

4.05 

165.62 

234.05 

4.12 

High-QS 

163.64 

248.90 

3.97 

165.44 

236.44 

4.03 

Combined data (0-150 GPa) 

163.81 

245.59 

4.02 

165.62 

233.93 

4.06 

Combined data (0-90 GPa) 

163.87 

243.41 

4.09 

165.67 

232.32 

4.11 

Combined data (0-75 GPa) 

163.93 

241.25 

4.15 

165.68 

232.29 

4.12 
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